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1 Introduction 



In this paper we shall explore Euclidean structures on manifolds which lead to 
Laplace operators. Euclidean structures can be introduced on a triangulation of 
a manifold by giving each simplex the geometric structure of a Euclidean sim- 
plex. This structure gives the manifold a length space structure in the same way 
a Riemannian metric gives a manifold a length structure: the length between 
two points is the infimum of the lengths of paths between the two points. The 
length of a path is determined by the fact that each simplex it passes through 
has the structure of Euclidean space. 

The purpose of this paper is to be able to do analysis on the piecewise Eu- 
clidean space. The Laplace operator A is well defined on many geometric spaces, 
and is especially important as a natural operator on a Riemannian manifold and 
as a generator of Brownian motion. In this paper, we define a general Euclidean 
structure called a duality triangulation which not only allows one to measure 
length between points and volume of simplices, but also allows one to describe 
a geometric dual cell decomposition and the volume of dual cells. This allows 
one to define a Laplace operator in a natural way, which has been applied to 
fields such as image processing HZ] and physics 36 . 

The duality triangulation structure is very similar to other Euclidean struc- 
tures used in both pure and applied math; specifically, we address the connection 
to weighted triangulations and Thurston triangulations. In addition, positivity 
of volumes of certain duals correspond to Delaunay or regular triangulations, 
which are used in a very wide range of applications from biology to physics to 
computer graphics. 

This paper is organized as follows. We begin in Section |21 with an intro- 
duction to Euclidean structures by recalling the definitions of weighted and 
Thurston triangulations, introducing dual triangulations, and relating the three 
types of triangulations. In Section |31 we discuss regular triangulations and De- 
launay triangulations and consider flip algorithms for constructing regular and 
Delaunay triangulations. In Section 0] we introduce the Laplace operator A 
associated to a given duality triangulation and derive some of its properties. Fi- 
nally, in Section El we briefly discuss the status of piecewise linear Riemannian 
geometry. 

The major new results in this paper are the result on the equivalence of 
weighted, Thurston, and duality triangulations in Section [2.51 the analysis of 
flip algorithms in Section |31 the generalization of Rippa's theorem to regular 
triangulations in Section ^^nd the definiteness results in Section 

Many of the results in this paper were motivated as generalizations of those 
described in 
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2 Euclidean structures 



2.1 Basic definitions 

In this section we shall introduce three types of Euclidean structures: weighted 
triangulations, Thurston triangulations, and duality triangulations. All struc- 
tures begin with a topological triangulation T = {Tq, 7i, . . . , 7„} of an n-dimensional 
manifold (we shall usually use n to denote the dimension of the complex in this 
paper). The triangulation consists of lists of simplices a'^, where the super- 
script denotes the dimension of the simplex, and 7^ is a list of all fc-dimensional 
simplices a'' — {iq, ■ • ■ , ik}- Wc shall often refer to 0-dimensional simplices as 
vertices, 1-dimensional simplices as edges, 2-dimensional simplices as faces or 
triangles, and 3-dimensional simplices as tetrahedra. We shall often denote ver- 
tices as j instead of {j} . Let denote the directed edges, where we distinguish 
from When the order does not matter, we use {i,j} to denote an 

edge. A triangulation is said to be an n-dimensional manifold if a neighborhood 
of every vertex is homeomorphic to a ball in R". A two-dimensional manifold 
is often referred to as a surface. Throughout this paper we will be dealing 
exclusively with triangulations of manifolds or parts of manifolds. 

In order to give the topological triangulation a geometric structure, each edge 
is assigned a length £ij such that for each simplex in the triangulation there 
exists a Euclidean simplex with those edge lengths. We call such an assignment 
a Euclidean triangulation (T, £) , where we think of ^ as a function 

i:Ti^ (0,oo). 

The conditions on £ include the triangle inequality, but there are further re- 
strictions in higher dimensions which ensure that the simplices can be realized 
as (non-degenerate) Euclidean simplices. The restrictions can be expressed in 
terms of the square of volume, which can be expressed as a polynomial in the 
squares of the edge lengths by the Cayley-Menger determinant formula. Each 
pair of simplices tr" and (T2 connected at a common boundary simplex ct"~^ 
is called a hinge. In a Euclidean triangulation every hinge can be embedded 
isometrically in M". 

Euclidean triangulations have the structure of a distance space with an in- 
trinsically defined distance. Given any curve 7 whose length can be computed 
on each Euclidean simplex, we can compute the total length of the curve L (7) 
as L (7) — J2a ^ ^) where (7 n a) is the length of the curve in the 
simplex a (if the curve intersects the simplex many times, we simply add the 
contributions of each piece of the intersection). In particular, we can consider 
curves which are differentiable when restricted to each simplex (these are called 
piecewise differentiable curves). The intrinsic distance is defined as 

d (P, Q) = inf {L (7) : 7 is a path from P to g} . (1) 

The class of paths can be either taken to be piecewise differentiable or piecewise 
linear since length is minimized on piecewise linear paths, as explained in |50[ 
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Figure 1: A weighted or Thurston triangulation with corresponding circles at 
the vertices. 



Section 2]. A path which locally minimizes length is called a geodesic and one 
which globally minimizes is called a minimizing geodesic. 

We are now ready to introduce more structures on Euclidean triangulations. 

2.2 Weighted triangulations 

We begin with weighted triangulations. 

Definition 1 A weighted triangulation is a Euclidean triangulation (T, £) to- 
gether with weights 

We think of the weight Wi as the square of the radius of a circle centered at 
the vertex i. These weighted triangulations are used in the literature on regular 
triangulations such as |15j and Thinking of the weights in this way, in 
each n-dimensional simplex there exists an {n — l)-dimensional sphere which is 
orthogonal to each of the spheres centered at the vertices (this means they are 
perpendicular if they intersect, or else orthogonal in the sense described in j40| 
Section 40]). In this way, each simplex a has a corresponding center C (cr) , 
which is the center of this sphere, and the center has a weight wc{a) which is 
the square of the radius of this sphere. See Figures ^ and [3 

An important particular case of weighted triangulations is that when Wi — 
for all vertices i. This is the basis for Delaunay triangulations, but may not 
satisfy the Delaunay condition. We shall revisit this in Section O 
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Figure 2: A weighted or Thurston triangulation with corresponding spheres at 
the vertices. 

2.3 Thurston triangulations 

Definition 2 A Thurston triangulation is a collection (T, w, c) , where 

w : To -> M, 
c : Ti ^ M, 

where Cij < Wi + Wj and such that the induced lengths 

make (T, t) into a Euclidean triangulation. 

For a Thurston triangulation, one considers the weight Wi to be the square of 
the radius of a sphere centered at vertex i, just as for weighted triangulations, 
and one considers Cy = 2rirj cos {tt — 9ij) where 9ij is the angle between the 
spheres centered at vertices i and j. In this case, one derives the formula for lij 
by the law of cosines. By considering Cy instead of Oij , we have included some 
cases where the spheres do not intersect. These structures were studied by W. 
Thurston in the context of proving Andreev's theorem (see |H2 and 

An important special case is that when Cy = —2rirj (i.e. Oij = 0). This is 
the case of a sphere packing on each simplex, since it corresponds to the spheres 
being mutually tangent (as in 1211] [21] )■ 
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2.4 Duality triangulations 

Definition 3 A duality triangulation is a collection (T, d) , where 

d : ^ R 

which satisfies 

dl + d% + dl = d% + dlk + dlj (2) 
for each {i, j, fc} e 72 and such that the induced lengths 

£{j — dij ~t~ ^ji 

make (T, i) into a Euclidean triangulation. 

We think of the weight dij as representing the portion of the length £ij of 
edge {i,j} which has been assigned to vertex i while dji is the portion assigned 
to vertex j. We thus call them local lengths. The total length of {i,j} is the sum 
of the contributions dij from vertex i and dji from vertex j. Hence each edge 
is assigned a center C {{i,j}) which is distance dij from vertex i and distance 
dji from vertex j. The condition |(2Jl ensures that for each triangle k} , the 
perpendiculars to the three edges through the edge centers meet at one point, 
which can be called the center of the triangle, C {{i,j, k}) . We shall soon see 
that this condition on 2-dimensional simplices allows us to define a center for 
every simplex in the triangulation. 

There are two canonical examples which automatically satisfy the condition 
(|2Jl. One is the case where dij depends only on i for all edges {i,j) (that is, 
dij = dik, etc.). We call this a circle or sphere packing as in j20) . and the dual 
comes from the inscripted circle, that is, the center C fc}) is the center 
of the circle inscribed in fc} in 2D and the center C ({i, j, A:,^}) is the 
center of the sphere tangent to each of the edges of the tetrahedron {i, j, fc, £} in 
3D. Another important case is where dij — dji. This corresponds to the center 
C ({«, j, fc}) coming from the circle circumscribed about the triangle {i, j, A;} and 
similar for all higher dimensions. 

The structure is called a duality triangulation because the existence of a 
center C {a) for each a puts a piecewise-Euclidcan length structure on the dual 
of the triangulation in such a way that dual simplices are orthogonal to ordinary 
simplices. For example, in two dimensions, if an edge is part of the 

two simplices fc} and then we can define the length of the dual 

edge ^{i:j} to be equal to the distance from the center C({z,j, A:}) of the 
triangle {«, j, fc} to the center C ({i, j}) of the edge {«, j} plus the distance from 
C({i,j,-^}) to C({z,j}). When the hinge is isometrically embedded in M^, we 
see that is a straight line which is perpendicular to the edge 

We shall now show that this can be done in all dimensions, and no additional 
restrictions must be made besides Q for each triangle. 

Proposition 4 A duality triangulation in any dimension has unique centers 
C (cr™) for each simplex cr™ such that C (cr™) is at the intersection of the 
(rn — 1)- dimensional hyperplanes through C{{i,j}) and perpendicular to {i,j} 
for each {i,j} in cr™. 
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Proof. We construct the centers C (ct"*) inductively for m-dimensional sim- 
plices. Each pair of m-dimensional simplices meeting at an (m — l)-dimensional 
simplex (a "hinge") can be embedded in as two adjacent Euclidean sim- 
plices. To make the notation more readable, we shall not distinguish between 
the embedding of the hinge in M™ and the hinge as abstract simplices in the 
piecewise Euclidean manifold. A simplex cr™ is assumed to be Euclidean with 
the assigned edge lengths given by . We now inductively construct the centers 
of each simplex. First, C {{i}) = i and C{{i,j}) is the point on which 
is a distance c?,y to {i} and a distance dji to {j} . Now, given centers C (cr*^) 
for A; < m — 1, we construct C {a™) as follows. Label the vertices of cr™ to be 
{0,1,..., m}. 

Let ^{ij} denote the plane in R™ through C and perpendicular to 

{i,j} (this is a hyperplane in M™). First we construct the center of a simplex 
{0,1,2} (m = 2). One can embed the simplex in as the three vertices 
(0, 0) , (^01, 0) , and {£02 cos 70, io2 sin 70) , where 70 is the angle at vertex 0. The 
centers of the three edges are realized as C({0,1}) = ((ioi,0), C({0,2}) = 
(do2 cos 70, do2 sin 70) , and C ({1, 2}) = {ioi - di2 COS71, di2 sin 71) . Hence 

n{o,i} = {(doi,i) :t GM}, 

n{o,2} = {(^02 COS70 -f tsin7o,do2sin7o -icos7o) : t G R} , 
n{i,2} = {(^01 - di2 cos 71 -I- i sin 71 , c?i2 sin 71 -f t COS71) : t g R} . 

A quick calculation (using the law of cosines to compute cos7i and sin7i in 
terms of dij ) shows that the three intersection points of these lines coincide if 
and only if (O holds. 

We now construct C (cr™) given C (a™~^) for all (m — l)-dimensional sim- 
phces. Since cr™ is a nondegenerate Euclidean simplex, the planes n|o,i}, • . ■ , ^{o,m} 
intersect at one point, c. We need only show that the planes ^{ij} also inter- 
sect c. This is true because inside {0,i,j}, the planes Hj-Q j}. and ^{o,j} meet 
each other and the plane ^{i^} at C ({0, i, j}) . Furthermore, since these planes 
are all perpendicular to {0,i,j} , the intersection Hj-Q jj n ^{ij} is equal to the 
intersection H^o.i} H H^o.j} and hence contains c. We call this point C (cr™) = c. 
■ 

Centers allow a geometric description of the Poincare dual of the triangula- 
tion. Any triangulation of a manifold has a cell complex which is its Poincare 
dual (see, for instance, [3| or HD). As noted by Hirani |23|, the assignment of 
a center to each simplex allows one to assign a geometric Poincare dual, or just 
dual for short. See Figures |31 and ^ for two-dimensional and three-dimensional 
simplices with dual cells included. Hirani restricted himself to "well-centered" 
triangulations, which means that the center of each simplex is inside the sim- 
plex. This is a very strong restriction, for even Delaunay triangulations may 
not be well-centered. Duality structures allow one to define geometric duals (a 
realization of the Poincare dual), each of which has a volume. The structure 
may not be well-centered, and for this reason some volumes may be negative. 
The fc-dimensional volume of a simplex cr'' will be denoted |(t''| (for instance 
|{*;i}| = ^ij) and the (n — fc)-dimcnsional (signed) volume of the dual of a 
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Figure 3: Two triangles with the pieces of dual edges intersecting the triangles 
included. 



simplex ifcr'^ will be denoted ■ 

It is helpful to consider an example before considering the general definitions. 
Given a triangulation of a three-dimensional manifold, one defines the duals as 
follows (compare with Figure^: 

0. The dual of a 3-simplex k, i} is the center, if k,£} = C ({«, j, k, £}) , 
and its volume is one. 

1. The dual of a 2-simplex k} contained in k, £} and {i, j, k, m} is 
a 1-cell if{i,j,k}, which is the union of the line from C{{i,j,k,£}) to 
C {{i,j,k}) and the fine from C {{i,j,k,m}) to C {{i,j,k}). Its volume is 
slightly tricky. We define the volume as 

I ★ {z, J, A;} I = ±d [C k,£}),C{{i, J, k})] ±d[C k,m}),C k})] 
= ±d[Ci{i,j,k,£}),C{{i,j,k,m})] 

where d is the Euclidean distance in R'^ (these are well defined because 
we can embed the hinge in R^) and the signs are defined appropriately. 
In the first fine, the sign is positive if C {{i,j,k,£}) is on the same side 
of the plane containing the side {i,j,k} as the simplex fc,£} is, and 
negative if it is on the other side (similarly for {i,j,k,m}). The sign on 
the second line is defined to be compatible with the previous definition. 
Note that it is possible for \ic {i,j,k}\ to be negative. 

2. The dual of a 1-simplex is the union of triangles. For each k,£ 
such that {i,j,k,£} is a simplex, the intersection of the simplex with the 
dual {i, j} is the union of the right triangle with vertices C {{i, j, fc, £}) , 
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Figure 4: Two tetrahedra with the pieces of dual edges and faces intersecting 
the tetrahedra included. 



C k}) , C and the right triangle with vertices C k, £}) , 

C i}) , C ■ Each of these triangles has a signed area. The first 
is 

±^d [C k, £}) , C {{z, J, k})] d [C ({z, j}) , C {{z, J, k})] 

and the second is defined similarly. The sign is defined as the product of 
the appropriate signs in each of the two distances. 

3. The dual of a vertex {i} is a union of right tetrahedra. For each j, k, I 
such that /c,£} is a simplex, the intersection of "A"!*} with 
is the union of six tetrahedra: 

(a) the tetrahedron defined by the vertices C ({«, j, fc, £}) , C({«,j, fc}), 
C({«:J}) : and i, 

(b) the tetrahedron defined by C ({i, j, fc, £}) , C({i,j, /c}), C({i,fc}), 
and i, 

(c) the tetrahedron defined by C ({i, j, fc, £}) , C I}) , C , and 

■i, 

(d) the tetrahedron defined by C k, £}) , C {{i, j, £}) ,C{{i, £}) , and 

(e) the tetrahedron defined by C {{i,j,k,£}) , C {{i,k,£}) , C{{i,k}), 
and «, 

(f ) and the tetrahedron defined by C {{ij, k, £}) , C ({i, k, £}) ,C{{i, £}) , 
and i. 
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The volume of {i} is the sum of the volumes of these tetrahedra, namely 
±^d [C {{t, J, k, £}) , C ({z, J, k})] d [C j}) , C k})] d [z, C 

for the first and similarly for the others, where the signs are defined ap- 
propriately. 

We can define the geometric duals in a triangulation of an rt-dimensional 
manifold inductively as follows. 

Definition 5 Define the dual of {0, ... , n} to be if {0, . . . , n} = C ({0, . . . , n}) , 
and |*{0,...,n}| = 1. 

Definition 6 The signed distance 

d± [CK),C(a"-i)] 

for (t"~^ C cr" is equal to the distance between C (fi") and C (cr"~^) in any 
isometric embedding cr" C K." with the sign positive if C {cr"') is on the same 
side of the hyperplane defined by u"^^ c K" as cr" is, and negative if C {a") is 
on the opposite side. 

It will be useful to know the following formula for the distance between the 
center of a triangle and the center of a side. Consider a triangle {«, j, fc} . Then 
some basic Euclidean geometry yields 

d± [C {{^,J, k}) , C ({z, J })] = (3) 

sm7i 

where 7^ is the angle at vertex i. 

Proposition 7 For any /c > 1, the volume of a simplex cr*' is 

^ n^±[^H'^(-^'^^)] (4) 

where is fixed and the sum is over all strings of simplices contained in . 

Proof. The proof is by induction on fc. If fc = 1, then |{i, j}| = dij + dji. 
Assume Q is true and consider a^^^ . Let the boundary of cr*^+-'^ be made up of 
ctq, . . . , cr^j^-^. The volume can be computed as 



i=0 



where each term in the sum is the volume of the simplex consisting of the 
center C (cr''+^) union cr*^ and the signs for d± tell us whether to add the area 
or subtract the area. It follows from the inductive hypothesis that 
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Note that the above argument works for any choice of center C (tr'^) G M*^ 
as long as C (cr^) are the orthogonal projections onto the subspaces spanned by 
cr^ for each subsimplex. The volume of a dual simplex is defined as follows. 

Definition 8 The volume of a dual simplex -j^a^ is defined to he 

n— 1 

= E \{d^[C{a^),C{a=^^)] (5) 

where is fixed and the sum is over all strings of simplices containing . 

Note that the volume is signed (it may be negative) . We note that the total 
volume is expressible in terms of volumes of the dual simplices. 

Proposition 9 Given a duality triangulation T of dimension n, the total vol- 
ume is 

Proof. We know that 

{i}C---Co-" j=0 

by © and 

-. n — 1 

by Q . Hence it is sufficient to show that 

E E 

ieTa {i}c---Ccr" 

is a reordering of 

E E • 

(T"er„ o-oc-'-cct" 

Here is one way to see this. Make a graph whose vertices are all simplices 
of all dimensions and whose edges connect two simplices if one simplex is in 
the boundary of the other. An easy way to draw the graph in the plane is to 
put vertices corresponding to n-dimensional simplices in a horizontal line on 
top, then (n — l)-dimensional simplices in a horizontal line below those, and 
so on until at the bottom is a horizontal line containing all of the vertices 
corresponding to 0-dimensional simplices in the triangulation. Now draw the 
edges, which can only connect a vertex in a row to a vertex in the row above 
or below. Now we shall see that both sums are equal to the sum over all paths 
between the top and bottom of this graph. We can count this in two ways, first 
start at the bottom with each path starting at a 0-dimensional simplex, or first 
start at the top with each path starting at an n-dimensional simplex. These are 
the two sums. ■ 
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2.5 Equivalence of metric triangulations 

We shall now show that weighted triangulations are equivalent to Thurston 
triangulations, and that, up to a universal scaling of the weights, both are 
almost equivalent to the set of duality triangulations. This is motivated by the 
geometric interpretations of the lengths, weights, angles, etc. 

First we show the equivalence of weighted triangulations and Thurston tri- 
angulations. 

Theorem 10 There is a bijection between weighted triangulations and Thurston 
triangulations. 

Proof. The definition of Thurston triangulation gives the map to weighted 
triangulations, keeping Wi the same and assigning 

iij = ^yw^ + Wj - Cij . 

Since we assumed that Wi + Wj — Cy > 0, £ij must be positive. Similarly, we 
can map the other way as 

Cij ^ Wi + Wj - 

Note that since £ij > 0, we must have that Wi + Wj — Cij > 0. ■ 

Next we map weighted triangulations to duality triangulations. Notice that 
there is a one parameter family of deformations of a given weighted triangulation 
of a triangle {i, j, k} which fix the center C {{i,j^ A:}). These deformations are 
given by 

Wi Wi+t (7) 
for varying t. We call these weight scaling deformations, or just weight scalings. 

Theorem 11 Weighted triangulations modulo weight scalings can be mapped 
injectively into the set of duality triangulations. It is a bijection if the set of 
duality triangulations are required to satisfy 

r 

V f d ■ , - ,, ) = (8) 

fc=0 

for all loops j = iq, ii, . . . , v = j , where {ik, ik+i} G Ti- 

Proof. The key observation is that given spheres at the vertices of a simplex 
with given radii ^ywi, one can always construct a sphere which is orthogonal 
to each of these spheres. The center of that sphere will be the center of the 
simplex, and for that reason is often called the orthogonal center 15 . By the 
arguments above, we need only construct the dual for triangles. One can do this 
very easily by embedding the circles in a vector space of signature 1,1,1,-1 as 
in |4U[ 40.2]. Given a center, one can draw the lines perpendicular to the sides of 
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the triangle through the center, and these determine dij. A careful calculation 
yields 

This is the map to duality triangulations. Note that the condition is auto- 
matically satisfied. 

There appears to be more information in weighted triangulations, however, 
because the new circle centered at the orthogonal center has a radius, which can 
be calculated to be 

dtk ~ dij cos7yfe^ ^ 



--'■'•■-^ r- —"'- I (10) 



ijk "'ij 

dfj + dfk - 2dijdik cos^ijk 



sin^ "/ijk 



Wi 



where jijk is the angle at vertex i in triangle {i, j, fc} . Note that rf^j. — wc({ij^k})j 
the weight assigned to the center of k} . The weight scalings allow, for any 
single triangle fc}, one to specify the value of rfjf, while fixing the center 
C k}) . Fixing the center means that each would map to the same duality 
triangulation. It is easy to see that the formula Q is unchanged by scaling de- 
formations like ((TJ. If one chooses r^-fc then the map is unique. Once this scale 
is fixed in one triangle, however, the scale is determined on adjacent triangles, 
because weights on shared vertices have been fixed, and the deformation |2J) 
must be done for all vertices i in the triangle. Thus there is one free scaling 
parameter for the whole triangulation (if it is connected). ■ 

The inverse map from duality triangulations to weighted triangulations must 
take dij + dji to £ij. In order to get the weights, we must first fix wq for a 
given vertex (this is a free parameter since we are considering the weighted 
triangulation modulo scaling) . Then each neighboring weight can be calculated 
using ©: 

Wj = d],- d^ +w.,. (11) 

We need only show that this is well defined. Suppose {i,j, fc} G 72 and consider 
a Wk which can be defined from Wj or Wi. Then we need that 

dl^ - dlk + = dij - rf^^ + Wj . 

But since Wj = dji — dfj + Wi, this follows from the fact that dl^ — d^^. = dij — 
d% + d% ~ dij from . It follows by a similar argument that any nuU-homotopic 
loop can be triangulated and property ||HJ) holds automatically, showing that for 
any nuU-homotopic loop j = iq, ?i, . . . , ?l = j of L vertices with {ik, ik+i} S Ti, 



y(dh ~d^ 



k=l 



Thus, in general, we need to assume property ((SJ is satisfied for the weights to 
be well-defined. For example, the following triangulation of the torus does not 
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Figure 5: A triangulation of the torus together with dual edges. 



satisfy ((HJ for all loops. Tile a torus with the two triangles {1, 2, 3} , {1, 2, 4} 
where d^i = ^21 = ^24 = 1 - e, rfis = c?i2 = ^42 = e, and c?32 = ^23 = 1^14 = 
dii = I for small £, see Figure El Note that 

dj^ + dl^ + = £2 ^ i + (1 - ef = dl^ + d% + 

d?2 + ^24 + dli=e'^ + \ + ef = dl, + dl^ + dl^ 

and so on. The homotopy-nontrivial loop containing {1,2} will not satisfy 
property ||SJ). However, if wc started with a weighted triangulation, property 
(|SJ) is automatically satisfied and thus the map from weighted triangulations to 
duality triangulations is injective. 

Corollary 12 For a triangulation of a simply connected manifold, there is a hi- 
jection between weighted triangulations up to scaling and duality triangulations. 

Proof. Since the manifold is simply connected, any loop bounds a 2-dimensional 

disk, homeomorphic to — G : \x\^ < 1^ , which is triangulated. One 

can easily prove by induction on the number of triangles triangulating the disk 
that on the boundary of any such disk, holds. ■ 

3 Regular triangulations 

3.1 Introduction to regular triangulations 

Recall the definition of a regular triangulation (see, for instance, ^Kj or 
Let d (x,p) be the Euclidean distance between points p and x. Define the power 
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distance 

TTp : M" M 

by 

np{x) ^ d{x,p) - Wp (12) 

if p is a point weighted with Wp. The power is important as a function which is 
zero on the sphere centered at p with radius y/w^, positive outside the sphere, 
and negative inside the sphere. Notice that if p is a vertex of a simplex a and 
c = C (a) then tTc (p) = Wp and vr^ (c) — Wc, where the weight Wc is defined as 
the square of the radius of the orthogonal sphere as described in Section 

Since we can embed any hinge in M", the following local definition of regu- 
larity makes sense on a piecewise Euclidean manifold. 

Definition 13 An{n~ 1) -dimensional simplex a"^^ incident on two n- dimensional 
simplices cr" ~ (t"^^U{vi} and — ct"^^U{u2} is locally regular ifuci {^2) > 
and TTcj (wi) > Wvi , where Ci = C (ct") is the center of cr" for i = 1 or 2. If 
the weights are all equal to zero, a locally regular simplex is said to be locally 
Delaunay. 

Sometimes we will instead say that the hinge is locally regular. A hinge is 
locally Delaunay if and only if it satisfies the local empty circumsphere prop- 
erty: the sphere circumscribing cr" does not contain V2- This is simply the 
interpretation of the definition when the weights are equal to zero. Note that 
the condition for being locally regular is unchanged by a weight scaling of the 
type Q due to the formula l(Tn|l for Wc({i.j,fc})- 

There are actually global definitions of regular and Delaunay, since the defini- 
tion of power H12|) makes sense globally using the intrinsic distance described 
m Section EH 

Definition 14 An n- dimensional weighted triangulation is regular if for every 
a" G 7^1, we have ttc{<t") (v) > Wy for every vertex v in the complement of a"". 
In the case that the weights are all zero, we say the triangulation is Delaunay. 

In the case of two-dimensional Delaunay, the condition on the power says 
that for every circle containing at least three vertices, there is no vertex inside 
that circle. It is a well known fact that for n-dimensional regular triangulations 
of points in R" [2] and for 2-dimensional piecewise Euclidean surfaces with 
zero weights [H] [HI] that every hinge being locally regular is equivalent to the 
triangulation being regular. It is likely that the proof in |32l Chapter 3] can be 
generalized to regular triangulations of any dimension, but we do not do that 
here. 

The argument in [2] uses the fact that a geodesic must be a straight line, 
and along a geodesic line the power increases in the manner listed below. To 
generalize that argument, one needs the following assumption: 

Criterion 15 Suppose the hinge {cr", j '''"^^j locally regular. Consider a 
minimizing geodesic ray 7 starting at Xq which intersects a hinge {cr", CT2 , ct"~^ } 
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by first entering a" and then The simplex u" ^ determines a plane which 
separates cr" and 02 and contains all points x such that '^(j(^„n'^ {x) — T^(j^„n'^ {x) . 

Then Trc{a^) (^0) < ^c{a^) (^0) • 

One might try to prove Criterion 1151 by "developing the geodesic" in the 
plane in the following way (we consider two dimensions for simplicity). Start 
with a triangle and embed it in M^. For each new triangle which the geodesic 
goes through, embed a copy in adjacent to the previous triangle so that it 
looks like we are unfolding the manifold. The geodesic must be a straight line 
if it does not go through a vertex and so we may try to make comparisons on 
this development. Note also that by the following theorem of Gluck, every two 
points have a minimizing geodesic between them. 

Theorem 16 Prop. 2.1]) If a piecewise Euclidean manifold is complete 

with respect to the intrinsic distance, in particular if M is a finite triangulation, 
then there is at least one minimizing geodesic between any two points of M. 

The problem with this is that geodesies do go through vertices and even by 
varying the endpoints slightly, a minimizing geodesic may still go through the 
vertex (see [SZl Figure 14]). Hence it is not at all clear that Criterion 1151 is 
always satisfied. 

Note that Bobenko and Springborn ^ are able to prove that Delaunay 
is the same as all edges being locally Delaunay in general by developing the 
triangulation (not along a geodesic). Their argument appears to strongly use 
the fact that the edges are locally Delaunay (with all weights equal to zero), 
but does not use Criterion El 

For completeness, we include the proof for regular triangulations of n-dimensional 
manifolds, assuming Criterion ll5l which is proven using a similar method. 

Theorem 17 Under the assumptiong of Criterion \15i an n-dimensional weighted 
triangulation is regular if and only if all of its hinges are locally regular. 

Proof. This proof is essentially the one seen in |2j for Delaunay triangulations. 
Clearly if the triangulation is regular, then all hinges are locally regular. Now 
suppose all of the hinges of a weighted triangulation are locally regular. Given 
a vertex v and a simplex cr" such that v is not in cr", we may consider the line 
L from ti to a point in the simplex cr". Possibly by adjusting the line slightly, it 
must intersect, in order, a sequence of n-dimensional simplices cr", . . . cr^? = cr" 
where w is in a simplex bordering cr". By Criterion ll5l we know that 

for « = l,...,fc— 1. Since the triangulation is locally regular, 

Wy < 7rp(^„) {v) . 
Stringing these together, we get that 

Wv < 7rc(<T") {v) ■ 
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Although we have not proven that regular triangulations and locally regular 
triangulations are the same, we will often suppress the word "local" in the rest 
of this paper, always considering the local property. 



3.2 Regular triangulations and duality structures 

In order to have a definition of locally regular in terms of duality structures, we 
first look at the two-dimensional case. A regular hinge {{i,j, k} , {i,j,i}} must 
satisfy 

7rc({i,j,fe}) W = d{C {{i, j, k}) , {£}f - rfj^ > m 
7rc({i,j,«}) {k) = d{C{{i,j,£}),{k}f - r?-^ > Wk- 

Proposition 18 The center C {{i,j,k}) and radius rijk are uniquely deter- 
mined by the three equations 

d (C k}) ,{i}f - rfj^, = Wi 

d{C{{i,j,k})dj}f-rlk=^3 
d{C{{i,j,k}),{k}f -rlk=Wk. 

Proof. Put the triangle in Euclidean space with vertices vi = 0,Vj,Vk- We 
know that C {{i,j, k}) = xvj + yvk for some x and y and let z be the unknown 
radius. Now we can write the first two equations as 

\xvj + yvkl"^ - z"^ = Wi 
\{xvj + yvk) - Vj\^ - = Wj 

so 

Wi - 2vj ■ {xVj +yvk)+ if J = Wj 
which is linear in x,y. Similarly, we have 

Wi - 2vk ■ {xVj + yvk) + ifk = Wk- 

So the problem reduces to a linear system 

Wi + 4 - Wj = 2£lx + 2£,j£,k (cos 7,) y 
Wi + ^ik -Wk = 2£ij£ik (cos 7i) X + 2£l^y, 

where 7i is the angle at vertex i, with solutions 

_ {wi + £lj - Wj) £ik - {wi + £f^ - Wk) £ij cos7i 
^~ 2 (sin^ 7,) il£,k 

_ {wi + - Wk) £ij - {wi + £lj - Wj) £ik cos7i 
^~ 2{sm^^i)£ij£l 

17 



and 



■ 

Corollary 19 IJ an edge is on the boundary o J regular, i.e. 

T^c[{i,j,k}) W =di.C {{i, j, k}) , - rjjf. = wi, 
then C k}) = C {{i,j,£}) and riju = riji. 

Proof. If d (C ({i, j, k}) , tf—rf^f^ = wi then (C ({«, j, fc}) , r^fe) satisfy the same 
three equations as (C £}) , which determine these uniquely. Hence 
they must be equal. ■ 

Corollary 20 An edge {i, j} is regular if and only if > 0. 

Proof. Clearly j}| = on the boundary of regular as in Corollary 1191 

since the centers are the same. It is clear that > if the edge is 

regular. ■ 

One can now address the case of n dimensions. The corresponding proofs 
go through essentially untouched, and one has the following characterization of 
regular triangulations. 

Proposition 21 An (n ~ 1)- dimensional simplex a'^~^ which forms a hinge 
with simplices tr" — a"^^ U {i} and a" = a"^^ U {j} is regular if and only if 
> 0. 

Note that •ka"~^ is a one-dimensional simplex, so the property of being 
regular has to do with lengths dual to {n — l)-simplices being positive. The 
previous discussion motivates the following definitions which, in light of Theo- 
rem ^2 S'l'G slight generalizations of those for weighted triangulations. 

Definition 22 An n-dimensional hinge at simplex a"^^ is said to be locally 
regular if |-^cr"^^| > 0. An n-dimensional duality triangulation T is said to be 
locally regular if > for all ct"^^ G 

The duality structure is called a Voronoi diagram in the case the triangula- 
tion is Delaunay. Voronoi diagrams can be described in a more direct way. A 
point x is in the Voronoi cell "A"!*} if it is closer to i than to any other vertex. 
The boundary of the Voronoi cells forms the (n — l)-dimensional complex called 
the Voronoi diagram. The analogue for regular triangulations is called a power 
diagram. A point x is in the power cell -A" {«} if its power distance tt^ [x) is less 
than TTj (x) for any j ^ i (see |2] El)- In the case of regular triangulations, the 
duality described in Section IT^ is the same as using power diagrams. However, 
our notion of duality is more general, making sense for weighted triangulations 
which are not regular. 

An interesting question is how to find a regular triangulation of a given 
manifold with given weights. One method of construction is via so called "flip 
algorithms." 
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3.3 Flips in 2D 

We first consider the case of two dimensions. One can imagine the following 
notion of a flip. Given a hinge consisting of two triangles {i,j,k} and 
incident on one common edge {i,j}, there exists a flip which exchanges this 
hinge with a new hinge, namely {i,k,£} and {j,k,£}. Note that the flip fixes 
the boundary quadrilateral which consists cyclically of the vertices i,k,j,£. This 
exchange is called a 2 ^ 2 bistellar flip, or Pachner move dSHl)- If the hinge is 
convex, then this can be done metrically. In fact, the flip can be made at the 
level of a duality structure. Given the hinge described above, to do the bistellar 
flip we need to construct dke and dgk such that the condition Q is satisfied 
in each of the new triangles. This is done by solving the following system of 
equations for dke and dik, 

d^k + dlt + d% = + d^, + 
dke + dik = d{k,£) 

where d (fc, £) is the distance between vertex k and vertex £. This distance is the 
Euclidean distance because the entire hinge can be embedded in M^. Note that 
the first equation is equivalent to 

d'jk + dli + djj = rffcj + + dJi^ 

using ||2Jl for triangles {«, j, k} and {i,j, £} ■ The system can actually be written 
in a form easier to solve: 

J I "fei + "if ^li "-ik /'1Q^ 
d{kj) ^^^> 

dke + dik = d{k,£) 

which is linear, although the dependence of d{k,£) on the remaining d's is not 
obvious (although easy to find using trigonometry). Hence the 2^2 bistellar 
flip is well defined on duality triangulations, and the triangle inequality follows 
automatically. The two hinges which are equivalent by bistellar flips are shown 
in Figure El 

The flip requires that the quadrilateral is convex, otherwise the flip would 
require that one part is folded back, which complicates matters. This motivates 
the following definition: 

Definition 23 A hinge is fiippable if the quadrilateral defined by the hinge 
when embedded in is convex. 

Now, given a convex quadrilateral, there exist two possible ways to make it 
into a hinge. The duals are uniquely determined by an assignment of centers to 
the edges on the quadrilateral. Let be the line perpendicular to {i,j} and 

through C ■ Then L^^ ^y and meet at a point which is the center 

C {{i, j, k}) and similarly L^^ gy and L^j iy meet at a point which is the center 
C{{i,j,£}). However, also L^i ^} and Ly^nj meet at a point which becomes 
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Figure 6: Two hinges differing by a bistellar flip, together with duals. 



C ({i, k,{\) after the flip, and similarly with L^j ^y and ^{j,^}- Hence the centers 
in the hinge form another quadrilateral dual to the hinge (see the right side of 
Figure EJ. One diagonal of the dual quadrilateral corresponds to and 
the other corresponds to -A" {fc, £} . One must have positive length and the other 
negative length (or both are zero if all dual lines meet at a single point), so either 
the hinge is regular, or it will become regular by a flip. One can also think of the 
flip of the hinge corresponding to a flip of the dual hinge. To make this argument 
rigorous, one simply uses the fact that {j, j} must be perpendicular to {i, j} , 
and considers the possible cases for being positive, negative, or zero. 

If it is negative, then it must look like the right side Figure Eland hence a flip 
makes |-^{fc,£}| positive. If is zero, then a flip maintains this. 

3.4 Flip algorithms 

The most naive flip algorithm is to take a given weighted triangulation, look for 
a fiippable edge which is not regular, and flip it. Continue until the triangula- 
tion is regular. This algorithm was first suggested by Lawson and shown to find 
Delaunay triangulations for points in (^3^, see also exposition in |2| and re- 
lated result in It was later shown to work for any 2D piecewise Euclidean 
triangulation (where the weights are all zero) independently in [2H1 and |47|. 
This turns out not to work to find higher dimensional Delaunay triangulations 
or to find regular triangulations (if there are nonzero weights) even in dimension 
2. It was later found that points in can be triangulated with regular trian- 
gulations (for any dimension) by incrementally adding one vertex at a time and 
doing all the flips before adding additional vertices. In this case one must pay 
close attention to the order of the flipping and the algorithm must either sort 
the hinges or dynamically decide which hinge to flip next QHI QH]' Unfortu- 
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nately, it is not yet clear how to extend these algorithms to piecewise Euclidean 
manifolds, since their proofs rely on the fact that the triangulations are in R". 
In this section we propose a subset of the space of all weighted triangulations 
for which the naive flip algorithm works, just as in the case of two-dimensional 
Delaunay triangulations. 

Consider the following set. 

Definition 24 A 2- dimensional duality triangulation is said to be edge positive 
if dij > for every directed edge (i,j) of the triangulation and for any possible 
flip, i.e. any solution of ^Iclp . 

Hence a triangulation is edge positive if the centers of each edge are inside 
the edge and if the center of the new edge after any flip is also inside that edge. 
This implies that any non-regular edge is flippable: 

Lemma 25 Given a 2D edge positive duality triangulation, if an edge is not 
regular, then it is flippable. 

Proof. We prove the contrapositive. Suppose a hinge consisting of {i,j,k} 
and {i,j,£} is not flippable, i.e. the quadrilateral is not convex. There can 
only be one interior angle larger than tt, and it must be at vertex i or j. Say 
it is at i. Let Lk be the line through vertex i which is perpendicular to k} 
and let Le be the line through vertex i which is perpendicular to {i,£}. Since 
dik > 0, the center C ({«, j, k}) must be on the side of Lk on which {i, k} lies; 
call this open half-space Hk. Similarly, C {{i,j,£}) must lie on the side of 
on which {i, £} lies; call this half space H^. Let Hj be the half-space containing 
{i,j} whose boundary is the line Lj perpendicular to {i,j} through i. Then 
C {{i,j, k}) must be in Hk n Hj and C {{i,j,£}) must be in Hk n Hg. Since Lk, 
Li, and Lj intersect at i and since the angle at i is more than tt, Hk H Hj and 
Hi n Hj are disjoint sectors in a half-space. Use Euclidean isometrics to make 
put the hinge such that i is at the origin, {i,j} is along the positive x-axis, and 
k has positive y- value (and hence £ must have negative y- value). Any possible 
segment must be on a vertical line which intersects {i,j}. It is easy 

to see that any such line must intersect Hk H Hj with a larger y-value than it 
intersects He n Hj, implying that {i, j}| > 0. ■ 

Theorem 26 The edge flip algorithm finds a regular triangulation given an 
edge positive duality triangulation. 

Proof. Since every flip maintains the edge positive property and every nonreg- 
ular edge is flippable, we can always do a flip if the triangulation is not regular. 
We now only need an monotone quantity which measures the progress of the 
algorithm to complete the proof in the same way as in ^S], |2H|i and j47| . 
Since we are in two dimensions, we can use the Dirichlet energy for almost any 
function, since the energy increases if a flip makes the hinge regular (see The- 
orem |^|. Since this function increases every time we perform a flip and there 
are finitely many possible configurations, the algorithm must terminate. ■ 
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Note that the edge flip algorithm to find Delaunay surfaces is a special case, 
since in that case, dij = £ij/2 > 0. In the next section, we suggest the analogue 
of this proof for higher dimensions. However, the analogue of edge positive is 
possibly less natural in this setting. 

3.5 Higher dimensional flips 

First let's consider the analogue of the 2 — > 2 bistellar move in higher dimensions. 
Recall that in any dimension, we can embed a hinge in K" . so the type of relevant 
flips must take place inside one or two simplices in M". The relevant flip is the 
2 ^ n flip in M" (see Figure H for the 3D version). The flip takes two simplices 
cr" = '^o"^ ^ {*} ^^"i "'7 ~ '^0^^ ^ {j} meeting at a common face ctq^^ = 
{fci, . . . , kn} and replaces it with n simplices a^^ = |i, j, /ci, . . . , fcp, . . . , , 

where kp indicates that kp is not present. The same argument as above shows 
that dij and dji can be chosen so that the duality conditions ((2Jl hold for each 
face and the choice is consistent because of the duality conditions which already 
hold. 

Now the duality structure gives a hinge a dual hinge similarly to above. Look 
at the Figure |S1 to see the 3D case. The boundary of cr" consists of the faces 

(T|5 = {fci, . . . , kn} and ct"^"^^ = |i, fci, . . . , fcp, . . . , /c„| for p = 1, . . . , n while 
the boundary of a" is similarly decomposed. Let L^n-i be the line through 
C ((t"~^) and perpendicular to ct"^^ for any {n — l)-dimensional simplex. We 
know that L^.^^ and La--^. intersect at the point C (cr^) for every p,q — 1, . . . , n 
by Proposition^ We can also consider after the 2 n flip. The boundary of cr^ 
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Figure 8: A flip in three dimensions together with dual ceUs. 



consists of cr"j, ^ and crjj, ^ together with cr^ j.^ = |z, j, A:i , . . . , fcp, . . . , fcg, . . . , fc„ | 

for q = 1, . . . , 71 and q ^ p. Hence L^n^^ and L^^^^ intersect at the point C (^<^kp^ 

for each p — 1, . . . ,n. We find that there is a polytope with vertices C (af ) , 

C (cr") , and C (^<^kp^ ^'^^ P — 1, ■ • ■ , This is the dual hinge. The centers 

C (cr") and C (c") are connected via the edge ★ctq^^ If |*cro"^| < then the 
flip on the hinge does a n ^ 2 flip on the dual hinge which results in removing 
ifaQ~^ and replaces it with "A'cr^'^,^ , which are (2) dual edges, each with positive 
length. 

We see that this sort of flipping is exactly what is needed to make regu- 
lar triangulations via some sort of flip algorithm. However, the condition of 
flippability is harder to guarantee. We now examine flippability. 

Definition 27 An n- dimensional triangulation is said to be m-central if C (cr*^) 
is inside cr^ for all k < m. 

So edge positive is the same as 1-central. Furthermore, n-central is what 
is called well-centered in |57]. We now show that (n — l)-central assures that 
nonregular hinges are flippable. 

Lemma 28 Given an {n — l)-central triangulation of an n-dimensional mani- 
fold, if a hinge is not regular, then it is flippable. 

Proof. The proof is essentially the same as the proof of Lemma Con- 
sider a hinge consisting of the simplices {z, fci, . . . , fc„} and {j, fci, . . . , A:„}. The 
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first claim is tiiat if the liinge is unflippable, tlien at least one dihedral angle 
must be greater than tt. This is clear because if every dihedral angle is less 
than or equal to tt, then the hinge is the intersection of half-spaces defined by 
the (n — l)-simplices on the boundary and hence convex. Now consider the 
hyperplanes whose dihedral angle is greater than tt. By relabeling we may as- 
sume that the hyperplanes are determined by faces (J^jT^ = {«, ^i, • • ■ , ^n-i} 
and crj^^ = {j, fci, . . . , and intersect at cto "^ = i^i' • • • ' ^n-i} • Because 

C (cr^"^) C cr^~^, the C (erf) must be inside the half-space defined by the plane 
Ilife^, the plane through ctq"^ and perpendicular to cr^~^, on the side containing 
(T^~^. We have the same for C (aj-) and since the angle is larger than tt we must 
have that [-A'ctq^^I > by a similar argument to that in the proof of Lemma 

EH ■ 

Regular triangulations of points in M" are usually produced via some sort of 
incremental algorithm (see JHli 12^1)' The key observation is that if a new point 
is inserted into a regular triangulation, then there is at least one non-regular 
hinge which is flippable (or there are no non- regular hinges and it is regular). 
The generalization to the manifold setting is the following. Let Star (v) , the 
star of a vertex w, be defined as all simplices containing v. 

Lemma 29 Suppose Criterion \15\ is true. If every hinge in a triangulation is 
regular except for hinges intersecting Star (v) for some vertex v, then some if 
some hinge is not regular, there exists a flippable nonregular hinge. Hence the 
triangulation can be made regular via a flipping algorithm. 

Proof (sketch). The proof in 15 (also with exposition in |14[ Section 12]) can 
be applied to this situation. We are able to prove this lemma in the generality 
of manifolds because we have supposed Criterion 1151 in that generality. ■ 

Using this lemma on subsets of R", one is able to construct regular triangu- 
lations by: insert one vertex, make the triangulation regular, and then insert the 
next vertex, make the triangulation regular, etc. Unfortunately, on a manifold, 
it is not clear what the intermediate triangulations are so the algorithm does 
not quite work. Also, if one starts with any triangulation, one may not have 
a regular triangulation which is reachable only by flips, as seen in the example 
Uni Fig. 5.1]. 

4 Laplacians 

Laplace operators on graphs and on piccewise Euclidean manifolds have been 
studied in many different contexts, for instance JOI) QOli 123 > pS]. 

PS) . |77) . 021 ■ The purpose of this section is to consider the comments 
from Bobenko and Springborn in [Sj, which suggests the use of Delaunay trian- 
gulations as a natural context in which to describe Laplace operators, and look 
at the generalization of these comments to regular triangulations. 
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4.1 Laplace operator defined 

The suggested Laplace operator on two-dimensional surfaces in [S] (also seen in 
|27|. \Sb\) is the following operator on functions f : Tq ^ M., 

(A/),;= J2 ih-.fi) (14) 

j:{jj}eri 

where Wij is defined by 

1 



y = o (cot7fc,j + cot -fuj) 



if 7feij is the angle at vertex k in triangle A:}, and the hinge containing 
consists of the triangles fc} and £}. Note that if Wij > then 
this is a Laplacian with weights on the graph defined by the one-skeleton of the 
triangulation, and that A/i > if fi is the minimal value of / and A/i < if /i 
is the maximal value of /. Bobenko and Springborn note that if the triangulation 
is Delaunay, then Wij > and the Laplacian is, in fact, a Laplacian on graphs 
in the classical sense (see |TT]'). 

A simple calculation shows that if we take the weights at all vertices to be 
zero, then the signed distance 

d± [C {{i,j,k}),C{{i,j})] = ryfeC0S7fe„ 

where rijk is the circumradius of triangle {i,j, k}. Since the circumradius can 
be computed to be 

2 sm 7fey 

we find that ^ 

d± [C {{i, j, k}) ,C{{i, j})] = -i^j cot j^j ■ 

It immediately follows that 

" " l{M-}l ■ 

We see that the Delaunay condition is equivalent to Wij > 0, which is equivalent 
to >0. 

In general, Hirani 1571 suggests the following definition of Laplacian: 

This formula has roots in the following integration by parts formula for the 
smooth Laplacian: 

A/ = / S/f -ndS (16) 

JdU 
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where n is the unit normal to dU. Taking 11 = ^ {i} and sUghtly rearranging 
terms, we get the corresponding formula on picccwisc Euclidean manifolds 



(A/), l*«l= E 



fj_ fi_ 

\{^,J}\ 



I*{^,J}I 



where 1 1^ ^. |'| is the normal derivative and {i, j}| is the surface area measure on 
the boundary of {i} ■ This formula is well defined on any duality triangulation 
(which is the motivation for the definition) and coincides with (|14|) in the case 
of Delaunay triangulations, except for the factor of One can think of 

the difference between considering the induced measure A/ dV instead of the 
pointwise Laplacian A/. It is, in fact, natural to consider the measure instead 
since, if we consider the discrete Laplacian approximating a smooth one, the 
pointwise Laplacian is only accurate when considered on scales larger than the 
scale of the discretization. 

We note that the Laplacian given by (|15|l is also the same as the Laplacian 
considered by Chow-Luo ID in two dimensions as observed by Z. He, where the 
duality is defined by Thurston triangulations as described above. It also appears 
in [20] [23 in three dimensions, where Thurston triangulations are considered 
such that dij depend only on i. Also, the Laplacian described in |33| is actually 
the Laplacian described above in (|14|l with the same weights Wij. The interest 
in these Laplacians is that they are not derived from means such as (llfill but 
instead as the induced time derivative of curvature quantities under geometric 
evolutions. 

The Laplacian defined in (|15|l is a Laplacian with weights on graphs in the 
usual sense (see m]) if the coefficients 



are each nonnegative. In two dimensions we see that this is implied by dij > 
and > 0, which is the condition that the triangulation is regular. 

Note that the Laplacian can be considered the gradient of a Dirichlet energy 
functional as described in 'H' , which is the analogue of the smooth functional 



The Dirichlet energy functional induced by the duality triangulation is 



This specializes in the case where the Wi — for all i £ Tq (or, equivalently, 
dij = dji = iij/2 for all € 71) to the Dirichlet energy in 0. Note that 

this energy is positive if j-A- {«, j}| > 0. 



l*«l 





(17) 
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4.2 A generalization of Rippa's theorem 

Rippa 0H| showed that if one considers the Dirichlet energy (|17|) on a triangula- 
tion of points in where the weights are zero (or equivalently, dy = dji — £ij/2 
for all edges flipping to make an edge Delaunay increases the Dirich- 

let energy. Bobenko and Springborn [S] note that his proof extends trivially 
to piecewise Euclidean surfaces (2-dimensional manifolds). We shall express 
Rippa's theorem in a way closer to the exposition on [B], which is in line with 
the notation in this paper. 

Theorem 30 ( | 46 | ) Let (T, £) be a piecewise Euclidean, triangulated surface 
with assigned edge lengths £, which we think of as a weighted triangulation with 
all weights equal to zero. Let Tq be the vertices of the triangulation and let 
f : Tq R be a function. Suppose T' is another triangulation which is gotten 
from T by a 2 2 bistellar flip on edge e (in particular, Tq — Tq,) such that 
the hinge is locally Delaunay after the flip. Then 

Et' if) < Et if) , 

where Eq- and Eq-i are the Dirichlet energies corresponding to T and T . As a 
consequence, the minimum is attained when all edges are Delaunay ( and hence 
the triangulation is a Delaunay triangulation). 

Rippa's proof involves calculating £'(/r') — £'(/r) and showing that it is 
negative. The key is a lemma which factors E (/rO ~ E ifq-) and for which we 
shall give a direct proof later for the more general case of regular triangulations. 
The only thing missing is the proof of the final sentence, which requires that 
flipping edges eventually produces a Delaunay triangulation, which is proved in 
|2HI and We can generalize the first part of Rippa's theorem to regular 

triangulations: 

Theorem 31 Let (T, d) be a duality triangulation of a surface with assigned 
local lengths d. Let Tq be the vertices of the triangulation and let f : Tq ^ 
be a function. Suppose (T',d') is another duality triangulation which is gotten 
from (T, d) by a 2 2 bistellar flip on edge e such that the hinge is locally 
regular after the flip. Then 

Et' if) < Er if) , 

where Eq- and Eq-' are the Dirichlet energies corresponding to (T, d) and (T', d') . 

The proof depends on the following important generalization of Rippa's key 
lemma ^BJ Lemma 2.2] (see also ^^). 

Lemma 32 Let T = {{1, 2, 3} , {1, 2, 4}} and T' = {{1, 3, 4} , {2, 3, 4}} be two 
hinges differing by a flip along {1,2}. Then 

E ifr') - E ifr) - ifr' (c) - fr ic))' Aj,,,^ 
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where 

2 [r^Ti - rir2) A1234 + Wi^234 + ^2^134 - ^3^124 - ^4^123 



8^123^134^234^124 



Aijk is the area of {i,j, k} , ^1234 = ^123 + ^124 = ^134 + ^234 is the area of 
the hinge, c is the intersection of the diagonals, ri is the distance between c and 
vertex i, and fT' and /r are the piecewise linear interpolations of f with respect 
to the different triangulations. One can write 

/r(c) = f^/2 + ^/i 

tl2 «12 

/t' (c) = ^/4 + Y^fz- 

«34 «34 

The proof is somewhat involved ahhough straightforward. We use a proof 
which is more direct than the ones given by Rippa 0^1 and Powar 0^ for the 
case of Delaunay triangulations. 

Proof. Because we are on a single hinge, it is equivalent to use weighted 
triangulations by Theorem 1111 Let {I, w) be the corresponding lengths and 
weights. A simple calculation tells us that 

d±iC {{i,j}) ,C {{i,j,k})) 1 , . , '^j . ^fc 
^ = 2 cot7fey + 2^ cot7j,fe + ^ cot7yfe - 



where ^ijk is the angle at vertex i in triangle {i,j, k} and Aijk = lii^ji k}\ is the 
area. For simplicity, we shall use the notation hij^k — d± (C {{i, j}) , C ({«, j, k})) , 
which we think of as the height of the triangle {i,j,C{{i,j,k})}. Note that 
\ic {1, 2}| = /ii2,3 + /ii2,4j for instance. For any function /, we can compute 



1 4 

Eifr')^Eifr) = -J2 «»^/»/^' 



2 



where 



ai2 = H — , ai3 



hl2,3 


^12,4 


1X2 


^12 


^14,2 


^14,3 


^14 


£14 


^24,1 


^24,3 


^24 


^24 



014 — —„ 7 , ^23 



O24 — -„ , 034 



"•13,2 "•13,4 

^23,1 ^23,4 

^34,1 ^34,2 

^34 44 



and an = — '^j^i o.ij (where we have symmetrized Oij = aji). We now wish to 
factor the coefficients. 

We can easily figure out in terms of areas in the following way. For a 
realization of the hinge, with Vi representing the coordinates of {i} , we see that 
c = ui + {v2 — vi) = V3 + (u4 — V3) . By taking the cross product with 
V2 — vi or V4 — V3 we find that 

^12^134 J ^34^123 

ri = — and r^ 



^1234 Ai 



234 
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where A1234 = A123 + Ai24 = Ais4 + A234 is the area of the entire hinge. 
Similarly, 

^12^234 J ^34^124 

r2 = — ; and n = —. . 

^1234 ^1234 

Thus 

fr, (c) - fr (c) = -/4 + ^/3 - ^/2 - ^/i 

= (^123/4 + ^124/3 — ^134/2 — ^234/1) • 

^1234 

Also useful will be the calculation 

r3r4 - rir2 = -j^ (^34^123-4i24 - ^?2^234^134) ■ 

^1234 

There are essentially two different types of coefficients to consider. We need 
only consider ai2 and 013 since the others are similar. Let ^ijk be the angle at 
vertex i in triangle k} . Consider ai2. 

/ll2,3 , /il2,4 

"''^^^^ 

1 Wi W2 W3 

= ^ cot 7312 + cot 7213 + 7775- cot 7l23 " 



2 '^'^ n\2 2ei2 4^123 

1 Wi W2 W4 

+ 2 + ^ W2 " 

1 / \ 1^1 / X 

= - (cot 7312 + cot 7412) + 77^5- (cot 7213 + COt72l4) 

Z 11-^2 

W2 , . . . \ ''^3 W4 



' 2£f,(^°*^-^ + ^°*'^-^) 4Ai23 Mi24- 

Let 6 be the angle at c in the triangle {1, 3, c} . We shall use the fact that in any 
triangle k} we have £ij = £ik cos jijk + ijk cos7j,/s to compute the parts. 

■^13^23 cos 7312 , ^14^24 cos 7412 
cot 7312 + cot 7412 = —. h 



2^123 2Ai 



24 



^Is --^12^13 cos 7123 ^14 - -^12^14 COS 7124 
2A123 2Ai24 

^13 , ^14 ff sin 7314 .a\ , f sin 7413 



- cotei + . , + cot6» 



2^4123 2^4124 \ \ sin 6* sin 7i23 / \ sin ^ sin 7124 

^13 _^ ^14 ^ f sin 7314 , sin 7413 



2A123 2^124 sin 6* \sin7123 sin 7124 

£h ^ ^14 1 42^134^1234 



2^123 2^124 sine 44^123^124 
^^3^124 + ^14^123 - ^^2^134 
2^123^124 

^13 + ^14 _j_ '^13--^124 + ^i4'il2:i ^ ' 12-4?34 ^12- Ulil- ^234 



2^1234 2Ai23^124^1234 2^123^124^1234 
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since 

sin 7314 = cos 7123 sin 9 + sin 7123 cos 9 

and 

sin 7413 = cos 7i24 sin 6 — sin 7124 cos 6. 

Furthermore, 

^?3^?24 + ^?4^?23 - ^?2^?34 = J^?2^?3^?4 (sm^ 7l24 + sin^ 7123 - sin^ (7123 + 7124)) 

= - ^^12^13^14 (sin 7123 sin 7124 cos 7134) 

= -2^123^124^3^14 cos 7134 



Since 

,2 A , „• 2 D „:„2 



sin"^ A + sin^ B - sin"^ (^4 + B) = -2 sin A sin B cos {A^-B). 
Thus we have 

(^13 + ^14 - 2^13^14 COS 7134) ^?2^134A234 



cot 7312 + cot 7412 = 



2^1234 2^123^124^1234 
34^123^124 — -^12^134^234 



2^1234^123^124 

^1234 



^123^124 

For the other parts, 



(r3r4 - rir2) . 



and 



Thus 



cos 7213 , COS 7214 
cot 7213 + cot 7214 = —. V —. 

sin72i3 sin 7214 
_ sin 7234 
sin 7213 sin 7214 

_ ^12^234 

2A123A124 



. , . ^12^134 

cot 7123 + cot 7124 = 



2A123A124 



1 W\ 

ai2 = - (cot 7312 + cot 7412) + (cot 7213 + cot 7214) 

+ ^(cot7l23 + COt7l24)-^-^. 



impUes that 



ai2 = -ri — ^/^^^/^^ . (2^1234 (r3r4 - rir2) + 1^1^234 + W2A134, - ti;3Ai24 - ^^4^123) 

4^123/1134^234^124 
= 2^234-4134$. 
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(Recall 

^ _ 2^1234 (?'3?'4 - rir2) + W1A234 + W2A134 - W3^124 - ^4^123 

8^123^134^234^124 

as in the statement of the lemma.) 
Now consider ais. We can compute 

/ll3,2 /ll3,4 



ai3 



«13 *13 

1 Wi W3 W2 

2 ^''' + Ws ''''' ^ 2^ ''''' " 41;^ 

1 Wi Ws 104 

- cot 7413 + -^TTT cot 7314 + TTTo- COt 7134 - 



134 



^ (cot 7213 - cot 7413) + (cot 7312 - COt 7314) 

W3 , V ^2 W4 



We see that 



2^,^(cot7l23-COt7l34)-4^^^^ , 

sin 7324 sin 7124 



cot 7213 - cot 7413 



Sin72i3smt' sin 7413 sin W 

^12^134^234 — ^34^123^124 



2^1234^123^134 

since sin 7324 = — cos 9 sin 7213+sin 9 cos 7213 and similarly sin 7124 = — cos 9 sin 7413+ 
sin cos 7413. We also get 

cos 7312 sin 7314 - cos 7314 sin 7312 



cot 7312 - cot 7314 = 



sm 7312 sin 7314 
sin 7324 



sin73i2 sin73i4 

^13^234 
2^123^134 



and 

, , ^i3-4i24 

cot 7123 - cot 7134 = — J— . 

ZA123A134 

And so 

— ^234^124 r A ^ \ a a a a \ 

ai3 = —. i—~A A — (2^1234 (r-3r4 - rir2) + W1A234 - ■W3A124 + t«2^i34 - ■W4A123) 

4^123^134 A234>ll24 
= -2^234^124$- 

A similar argument gives the other coefficients. Then we see, for instance, 
that 

ail = — ai2 — ai3 — ^14 

= 2 (-^234^134 + ^234^124 + ^234^123) $ 
= 2A^34* 
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with similar expressions for 022, 033, and 044. Finally, we get that 
E (/rO - E (/r) = (A123/4 + A124/3 - ^234/1 - A134/2)' 

which is equivalent to the lemma. ■ 

Now we can prove the theorem. 
Proof of Theorem I31L Note that since we are only concerned with a hinge, 
it is equivalent to consider weighted triangulations or duality triangulations. 
Since the coefficient 012 ~ ^*{i2}\^ ^^"^ ^34 ^ ^nTs^TT' ^ ^ 

and a34 < if and only if T' is regular after the flip and not regular before the 
flip. Since all areas Aijk are positive, ai2 < if and only if $ < and hence 
the result is proven. ■ 

Note that in the proof we have shown that <!> < if and only if T is not 
regular and T' is regular. 

In order to get the global statement, one needs to know that a regular 
triangulation can be found using flips. This is not true in general (see jl5)\ 
However, we investigated some conditions when a flip algorithm does work in 
Section in 

As a corollary of Rippa's theorem, we get an entropy quantity that increases 
under the action of flipping to make a hinge regular. 

Corollary 33 Consider the entropy defined by 

A = inf I : ^ /f = 1 and ^ = ol . 

Then A decreases when an edge is flipped to make the hinge regular. 

Proof. Let A' denote the entropy after the flip and let /o be the / which realize 
A (since / is in a compact set, there must be an actual / which minimizes E (/)). 
Then 

A' = inf £;t' (/) < Et' (/o) < Er (/o) = A. 

■ 

Note that A can be considered an eigenvalue of a particular operator closely 
related to A. We remark that Corollary is similar in spirit to what is proven 
by G. Perelman at the beginning of his paper (4I|, where he shows that a slightly 
more complicated entropy, 

infjy (^Rf+i\\/f\')dV: J = 

where R is the scalar curvature, increases under Ricci flow. 

Note that in n dimensions, the regularity condition corresponds to |-^cr"~^ | > 
while good Dirichlet energy corresponds to |'A'0'^| > 0. Hence the correspon- 
dence between regular triangulations and the Dirichlet energy only occurs in 
dimension 2 because 1 = 2 — 1, which is why the theorem is only described 
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for dimension 2. Although wc do not pursue it here, this may indicate that 
the Laplacian should instead be defined on functions on vertices of the dual 
complex, / : — > M, in which case the Laplacian would be 



(Af^ - ^ V 1^"^^°' (f -f ) 



where the sum is over all n-simplices. In this case, positivity of the coefficients 
corresponds to being regular. 

4.3 Laplace and heat equations 

Given a Laplace operator, we can now consider the standard elliptic and parabolic 
equations, namely the Laplace equation 

Am = (18) 

and the heat equation 

|.A„, (19) 

where the heat equation is an ordinary differential equation since A is a dif- 
ference operator. A solution u to the Laplace equation is called a harmonic 
function. 

In order to study these equations, it will sometimes be easier to consider 
Au = as a matrix equation. We think of u : To — > M as a vector and A 
corresponds to a matrix L whose off-diagonal pieces are 

and whose diagonal pieces are 



La-- ^ 



Then one can write the Laplace equation as 

Lu = 0. 

Note that if we wish to consider Poisson's equation 

Au = f (20) 

then this is equivalent to 

Lu = fV 

where (fV)^ = fi\ii:{i}\- It is clear that L has the constant functions fi = a 
(or the vector (a, a, ... , a)) in the nullspace. If {hj}\ > then we find the 
following. 
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Theorem 34 If\'fc {i, j}\ > for all edges {i, j} then L is negative semidefinite 
with nullspace spanned by the constant vectors. 



Proof. In this case we have an TV x TV matrix L with diagonal entries negative 
and off-diagonal entries positive and with X^jLi ^ij = 0- We reiterate an argu- 
ment from ^21- Let (wi, . . . ,vn) be an eigenvector corresponding to A > 0. We 
may assume that > is the maximum of Vi . We wish to show that Vi — Vj 
for all i,j. Observe 



Equality holds if and only if Vi vi for all i. ■ 

Corollary 35 // > for all edges {i,j} then Poisson's equation has 

a solution for any f such that 



This is the analogue of the smooth result that Am — f has a solution if 

fdV — 0. One may also consider boundary conditions such as Dirichlet and 
Neumann conditions. These cases for Delaunay triangulations in two dimensions 
were studied by Bobenko and Springborn (Hj. 

The condition > is obviously very important for the proof of 

Theorem 1341 In two dimensions, this condition is equivalent to being regular by 
Corollary|2ni It is not always necessary to assume {i, j}| > 0, as seen in the 
following special cases. 

Recall that in two dimensions, if a duality triangulation is edge-positive, then 
the flip algorithm finds a regular triangulation fThcorcnil26(). For a similar set 
of two-dimensional triangulations, the Laplacian is negative semidefinite. 

Theorem 36 For any triangulation such that dij > for all {i,j) G , the 
Laplacian matrix L is negative semidefinite with nullspace spanned by the con- 
stant vectors. 

We begin with a series of claims and an important lemma before begin- 
ning the proof We shall prove this by a sequence of claims. For all of the 
claims it is assumed that the weights dij are all positive. We shall use hij = 
d± [C({1,2,3}) ,C({i,j})] and is the angle at vertex i. Consider only the 
3x3 matrix M corresponding to {1,2,3} with entries Mij = hij/£ij \i i ^ f 
andM,, = ''Y.J^^^'hr 

Claim 37 If hij < then 7i < § and jj < ^. 

Proof. Let k be the third vertex so that {i,j, k} = {1, 2, 3} . We know that 



N 



N 




h,j = 



dik - dij cos 7j 



sin7i 
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by formula El . If hij < then < dik < dij cos Hence cos 7,; > and 
7i < tt/2. We can also express hij as 

_ djk - dji cos7j 
sm 7j 

and follow the same logic. ■ 

Thus only one Mij may be negative. Suppose it is M12. 

Claim 38 M12 + A/13 = ^^^(^12 00372+^13 cos 73) _ 
Proof. We calculate 

, (i23 -0^21 cos 72 d32-fi3lC0S73 

iWl2 + M13 = ; 1 : 

£i2Sm72 I;i3sm73 

_ ^23 (^23 - d2l COS 72 - dsi COS 73) 

~ 2A123 

and finally we use that ^23 = ^12 cos 72 + ^13 cos 72 . ■ 
Claim 39 di2 cos 72 + di3 cos 73 > 0. 

Proof. If both 72 and 73 are less than or equal to 7r/2 then this is clear (since 
both may not be equal to 7r/2). Since M12 < 0, and hence /112 < 0, we can only 
have 73 > tt/2. Since hi2 < and his > we have that 

di3 ^ ^ di2 

— < cos 71 < — 

«12 «13 

so di2 > c?i3. Furthermore, since 71 + 72 < tt we have that 
< — cos 73 = cos (71 + 72) < cos 72 

so 

-dl3C0S73 < di2COS72. 



Lemma 40 Ma < 0. 

Proof. By the above argument, we know that Afn = — Afi2 — Afi3 < 0. Similar 
arguments hold for the other coefficients. ■ 

Proof of Theorem I36L It is sufficient to prove that for any matrix My, 
1 < j < 3, is negative semidefinite. We know that the vector (1, 1, 1) is in the 
nullspace and we have already shown in Lemma 1401 that the diagonal entries 
are negative. Hence it is sufhcient to show that the determinant of the 2x2 
submatrix Mij, 1 < «,j < 2, is positive. We find that the 2x2 determinant 
is equal to Afi2Afi3 + Afi2M23 + A/i3Af23- We compute the determinant to be 
equal to 

(^13^23 + £^23^13) sin 72 

^12^3 
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(to do this calculation, begin by writing the terms in the determinant using 
formula JJJ choosing all of the denominators to contain sin 71 sin 72 , then re- 
arrange the terms using the facts that 71+72+73 — tt, dij + dji — £ij, and 
^ij = £ik cos ji+ijk cos7fe several times and finally recollecting /123 and his again 
using formula l|3Jl)- Note that the determinant is symmetric in all permutations 
in 1, 2, 3. We know by the claim above that two of the three hij must be posi- 
tive, so choosing the two that are positive, we must have that the determinant 
is positive. Hence the matrix is negative semidefinite. ■ 

We consider dij to be the length of a vector located at i and in the direction 
towards j. Thus the condition dij > is like a positivity (or Riemannian) 
condition for a metric (which measures the length of vectors) and is thus a 
somewhat natural condition. The following is another result on definiteness of 
the Laplacian with different assumptions. 

Theorem 41 For a three-dimensional sphere packing triangulation, L is nega- 
tive semidefinite with nullspace spanned by the constant vectors. 

Proof. It is proven in 2] (see also jlH]) that the matrix A/i 2,3,4} = ( ) 

is negative semidefinite with nullspace spanned by the vector (ri, . . . , r4). If we 
let i?{i,2,3,4} be the diagonal matrix with r^, i — 1, ... ,4 on the diagonal, we 
see that 

where (Mg.3)^ is the matrix extended by zeroes to a |7^| x |7^| matrix so that the 
(Afg-s)^ acts on a vector (ui, . . . , v\'Xo\) only on the coordinates corresponding to 
vertices in . Since > for all i e 7^, it follows that L is negative semidefinite 
with nullspace spanned by (1, . . . , 1) . ■ 

The importance of this result is it does not assume any positivity of the 
dual area, which appears to be stronger than the assumption that L is negative 
definite. If L is negative semi-definite with nullspace spanned by the constant 
vector (1, . . . , 1) then one can always solve the Poisson equation for / such that 

EM. = 0. 

The heat equation is an time-dependent, linear ordinary differential equation 

du 

— = Lu 
dt 

whose short time existence is guaranteed by the existence theorem for ordinary 
differential equations. One of the key properties of the heat equation is the 
maximum principle, which says that the maximum decreases and the minimum 
increases. This is true if j-A'lijj}! > 0. 

Theorem 42 // > then for a solution Ui (t) of the heat equation, 

Umax (i) decreases and Wmin (i) increases, where Umax = max{uj : i £ %} and 
Umin = min{ui : i e %} . 
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Proof. The proof is standard and is simply that for any operator Eu defined 

by 

for some weights > 0, then {Eu),- < if Ui = Mmax and {Eu)^ > if 

Note that the maximum principle is not equivalent to L being negative 
semidefinite; it is a stronger condition and the proof uses that the coefficients off 
the diagonal are positive. However, for certain functions (geometric ones which 
are related to the coefficients of the Laplacian), it may be possible to show that 
the maximum decreases and the minimum increases. We call this a maximum 
principle for the function / and we say that the operator is parabolic-like for 
the function /. In j^J it is proven that the sphere-packing case is parabolic-like 
for a curvature function K. 

5 Toward discrete Riemannian manifolds 

Much of this work arose out of an attempt to describe Riemannian manifolds 
using piecewise Euclidean methods. In this final section, we try to describe some 
of the work already done toward this end. There are two different philosophies. 
One is to find analogues of the Riemannian setting. The idea is to set up a 
framework on which variational-type arguments may be made analogously to 
those in the smooth setting. The other is to actually approximate smooth Rie- 
mannian geometry with discrete geometric structures. We shall briefly consider 
both of these. 

5.1 Analogues of Riemannian geometry 

In this paper we gave a discrete operator on duality triangulations which, it was 
argued, is an analogue of the Laplacian on a Riemannian manifold. This gives 
rise also to a discrete heat equation, which is an ordinary differential equation in 
this setting. It is not hard to imagine that similar arguments give rise to Laplace- 
Beltrami operators on forms with the proper deflnition of forms. A fc-form can 
be defined to be an element of the dual space to the vector space spanned by 
the fc-dimensional simplices. There are also dual fc-forms which are elements 
of the dual space to the vector space spanned by the duals of the {n — k)- 
dimensional simplices. Hirani |27j describes how to use duality information as 
we have described to define the Hodge star operation, and thereby the Laplace- 
Beltrami operator on these forms. One may then ask about an analogue of the 
Hodge theorem. This has been studied somewhat by R. Hiptmair |j26j. Study 
of the Laplace-Beltrami operator on manifolds is also related to the study of 
the Laplacian and harmonic analysis on metrized graphs and electrical networks 
(see [El, 0, H). 

Another important aspect of Riemannian geometry is the study of geodesies, 
which we recall are locally length-minimizing curves. In the setting of piecewise 
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Euclidean manifolds, the geodesies are piecewise linear. One may then ask 
many questions about geodesies, sueh as the number of elosed geodesies (see 
Pogorelov's work on quasi-geodesics on convex surfaces 021) and the size of 
the cut locus to a basepoint, the locus of points with two or more geodesies 
connecting it to the basepoint (see Miller-Pak ^7\). Many results on geodesies 
on piecewise Euclidean manifolds were found by D. Stone j5U) . which lead him 
to some possible definitions of curvature. The discrete geodesic problem for 
polytopes in R'^ was studied extensively in |38| . 

Much of modern Riemannian geometry is concerned with different notions 
of curvature, such as sectional, Ricci, and scalar. In the piecewise Euclidean 
setting, there are a number of definitions of curvatures, although it is still 
somewhat an open question which ones are the proper ones for classification 
purposes. Since the literature in this area is vast, we simply indicate some 
of the principle works. D. Stone jnO] was successful in proving analogues of 
the Cartan-Hadamard theorem (that negatively curved manifolds have univer- 
sal cover homeomorphic to M") and Myer's theorem (that positively curved 
manifolds are compact with a bound on the diameter) on piecewise Euclidean 
manifolds using a quantity which he calls bounds on sectional curvature. T. 
Regge introduced a notion of scalar curvature which is described at each (n — 2)- 
dimensional simplex as 27r minus the sum of the dihedral angles at that simplex 
This has been widely studied as the so-called "Regge calculus" (see, for 
instance, ^7], j^Sj, QUi DO)- There are even some convergence results, which 
we mention in the next section. Another potential curvature quantity in three 
dimensions is described by Cooper and Rivin in ^2). They consider the cur- 
vature at a vertex to be An minus the sum of the solid (or trihedral) angles at 
the vertex. This curvature is certainly weaker than the curvature introduced by 
Regge, but may be related to scalar curvature. It is possible that the right cur- 
vature quantity will lead to a geometric flow which simplifies geometry in a way 
similar to the way Ricci or Yamabe flow do in the smooth category. This has 
been studies a bit in ^HI) inHji 1201 > [211 1 and actually was the initial motivation 
for the definitions of Laplacian described in this paper. Other applications of 
discrete analogues of Riemannian geometry or geometric operators can be found 
in ini, |2E|) innii 021 ^ and [S2I. In addition, techniques applying to metric 
spaces with sectional curvature bounded in the sense of Alexandrov may apply 
(see i). 

5.2 Approximating Riemannian geometry 

Another goal is to approximate Riemannian geometry by a discrete geometry 
such as piecewise Euclidean triangulations. One would hope to be able to find 
elements of Riemannian geometry such as Laplacian, Levi-Civita connection, 
sectional curvature, scalar curvature, and so forth and not only have analogous 
structures, but be able to show that as the triangulation gets finer and finer, 
the discrete versions converge to the smooth versions. We mention here some 
of the results which have been successful in this direction. 

One of the most influential works is by Cheeger, Miiller, and Schrader, who 
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were able to relate discrete curvatures to Lipschitz-Killing curvatures [H]. The 
relevant discrete curvature is the sum certain angles and volumes of hinges. 
In particular, the scalar curvature measure {RdV) is concentrated on (n — 2)- 
dimensional hinges in a triangulation, and under a condition that the triangu- 
lation does not degenerate, they find that the curvature quantity 27r minus the 
sum of the dihedral angles multiplied by the volume of the {n — 2)-dimensional 
hinge converges to the scalar curvature measure. This version of scalar curva- 
ture is also the one suggested by Regge [321 and used extensively in the Regge 
calculus. They prove convergence for each of the Lipschitz-Killing curvatures. 
In addition, Barrett and Parker |S] proved a pointwise convergence of piecewise- 
linear approximations of the Riemannian metric tensor and certain types of 
tensor fields. 

In regards to the Laplacian, some experimental work has been done by G. 
Xu studying pointwise convergence of different discretized Laplace-Beltrami op- 
erators to the smooth ones [^^1 Some of the discretizations are the same 
or similar to those considered in this paper, while some are not. On graphs 
(one-dimensional manifolds and generalizations), it has been shown that the 
eigenvalues of the discrete Laplacians on metrized graphs converge to the eigen- 
values of the smooth Laplacian on a metrized graph JS] 

It was W. Thurston's idea to approximate the Riemann mapping between 
subsets of C by mappings of circle packings. Such a discretization has been 
shown to actually converge to the Riemann mapping |49| . 

Acknowledgement 43 / would like to thank Herbert Edelsbrunner, Feng Luo, 
and Igor Pak for some very helpful conversations related to parts of this paper. 
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